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Let  X(t),t» 0,  be  a  stationary  Gaussian  process,  and  define  the  sojourn  time  L„(r)  = 
mesfs:  0sss(,X(s)>u}  and  the  maximum  Z(t)  =  max(A'(s): 0«s  « /).  Limit  theorems  for  the  distri¬ 
butions  of  L„(f)  and  Z(t),  for  t,  u-»oo,  are  obtained  under  specified  conditions  on  the  spectral  density 
of  the  process.  The  results  supplement  earlier  theorems  obtained  under  suitable  conditions  on  the 
covariance  function. 

AMS  1980  Subject  Classifications:  60F05,  60G10,  60G15. 

stationary  Gaussian  processes  *  spectral  density  function  *  mixing  condition  *  sojourn  above  a  ievei  * 
extreme  value 


1.  Introduction  and  summary 

Let  X(t),  (?0,  be  a  stationary  Gaussian  process  with  mean  0,  variance  1  and 
covariance  function  r(t)  =  EX(0)X(t).  There  is  now  an  extensive  literature  on  the 
functionals  of  this  process  which  are  associated  with  the  large  values  of  the  sample 
functions.  Key  examples  of  such  functionals  are  (i)  the  sojourn  time  above  a  high 
level  u. 


Lu(t)  —  l[x(s)>ii]  ds; 

Jo 

. 

(1-D 

and  (ii)  the  maximum  value 

Z(t)  -  max  X(s). 

(1.2) 

0 


The  latter  is  well  defined  because  the  conditions  assumed  in  the  hypotheses  of  the 
theorems  imply  the  continuity  of  the  sample  functions.  Results  of  much  mathematical 
interest  are  limit  theorems  for  the  distributions  of  these  functionals  for  t  -*  oo,  u  -*  <x>, 
where  t  and  u  are  tied  together  by  means  of  a  specified  asymptotic  relation.  Among 
the  recent  results  for  Lu(t)  are  those  of  Berman  (1980,  1989).  Surveys  of  results  for 
Z(t)  are  those  of  Leadbetter  et  al.  (1983)  and  Leadbetter  and  Rootzen  (1988). 

This  paper  represents  results  obtained  at  the  Courant  Institute  of  Mathematical  Sciences,  New  York 
University,  under  the  sponsorship  of  the  National  Science  Foundation,  Grant  DMS  8801188,  and  the 
U.S.  Army  Research  Office,  Contract  DAAL  03-89-K-0125. 
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Central  Limit  Theorems  for  Extreme  Sojourns  of 
Stationary  Gaussian  Processes 

SIMEON  M.  BERMAN 
Courant  Institute 

This  paper  is  dedicated  to  Natascha  A.  Brunswick  for  years  of  talented  service, 

dedication,  and  good  cheer. 

Abstract 

Let  X(t),  t  s  o,  be  a  real  stationary  Gaussian  process,  and,  for  u  >  0  and  t  >  0,  let  L,(u)  be  the 
time  spent  by  -Y(r),  0  £  s  £  t,  above  the  level  u.  Here  u  is  taken  to  be  a  function  u(t)  of  t,  and  L,  is 
defined  as  L,(u(t)).  It  is  shown  that  the  distribution  of  (L,  -  £L,)/(Var  L,)l/:  converges,  for  r  — •  oo 
and  u(t)  —  co,  to  a  standard  normal  distribution  under  various  conditions  relating  the  growth  of  u(t) 
to  the  decay  of  the  covariance  and  other  functions  associated  with  it. 

"'3S?Sv  •  . 

Introduction  and  Summary 

Let  X(t),  t  ^  0,  be  a  real  measurable  stationary  Gaussian  process  with  mean  0 
and  covariance  function  r(t)  =  EX(0)X(t).  For  simplicity  we  take  r(0)  =  1.  For 
/  >  0,  put  L,(u)  =  mes(s  :  0  £  s  £  t,  X(s)  >  u);  then,  for  a  given  measurable 
function  u(t),  we  define 

(  Ll)  L,  =  L,(u(t))  -  f  1  [,r(j)>u(t;i  ds . 

Jo 

The  main  results  of  this  paper  are  several  new  central  limit  theorems  for  the  dis¬ 
tribution  of  (I,  -  ££,)/(  Var  L,)'12,  for  /  -»»  oo,  under  various  conditions  on  the 
covariance  function  and  the  function  «(f).  In  this  work  we  assume  that  the  spectral 
distribution  function  in  the  representation  of  r(t)  is  absolutely  continuous.  Let/(X) 
be  the  spectral  density  function,  and  let  b(t)  be  the  Fourier  transform  of  the  L2- 
function  (/(X))1/2;  then  r(t)  has  the  representation 

(1.2)  r(t)  =  |  b(t  +  s)b(s)  ds , 

and  b  E  L2.  Furthermore,  X(t)  has  the  stochastic  integral  representation 

( 1-3)  X(t)=  r  b(t  +  s)£(ds) , 

J— oo 

where  £(s)  is  a  standard  Brownian  motion. 

The  focal  point  of  the  calculations  in  the  proofs  of  these  theorems  is  the  relation 
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(2/m)  1/2  -*•  0,  which  proves  (6.7).  In  [2]  conditions  (1.5)  and  (1.7)  above  are 
used  to  prove  that  v(t)  -*•  oo  and  exp{5u2}/?(u/2)  -*•  0.  We  will  show  that  the 
latter  relations  remain  true  under  (1.5)  and  (1.7)  with  the  new  definition  (7.1 )  of 
v.  Indeed,  if  ( 1 .5 )  holds  for  some  9  >  1,  then,  for  every  9',  with  1  <  6'  <  6,  by  the 
relation  (5.6)  of  [2]  with  9'  in  the  place  of  9, 

v  ^  constant(t/«2)1/2exp{—  u29'/4} 


=  constant  tl/2exp{-«20/4} 


exp{(M2/4)(fl-fl')} 

u 


» 

; 

i 


from  which  it  follows  by  (1.5)  that  v  -*•  oo  and  by  ( 1.7)  that  e.xp{5u2}(3(v/2)  -*■ 
0.  In  the  relation  (5.9)  of  [2]  the  factor  mi/4  should  be  changed  to  u~U2. 

I  thank  Katarzina  Pietruska-Paluba  for  noticing  this  error. 
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A  CENTRAL  LIMIT  THEOREM  FOR  INTEGRAL  FUNCTIONALS  OF  A 
STATIONARY  GAUSSIAN  PROCESS* 

Simeon  M.  Berman 

Courant  Institute  of  Mathematical  Sciences 
New  York  University- 
New  York,  NY  10012 

ABSTRACT.  Let  X(t),  t  >  0,  be  a  real  stationary  Gaussian  process  with  covariance 
function  r(t).  Let  f(x)  be  a  function  in  £2(^)5  where  <j>(z)  is  the  standard  normal  density, 
and  assume  that  f  x  f(x)  (f>(x)  dx  0.  It  is  shown  that  the  central  limit  theorem  holds  for 
the  functional  fQ  f(X(s))  ds,  for  t  —*  oo,  under  the  sole  assumptions  r(t)  >  0  and  r(i)  — *■  0 
for  t  — ►  oo. 

1.  Summary. 

Let  X(t),  t  >  0,  be  a  real  stationary  Gaussian  process  with  mean  0  and  covariance 
function  r(t)  =  EX(0)  X(t).  For  a  Borel  function  f(x)  and  t  >  0,  consider  the  functional 

(1.1)  j  f(X(s))  ds  . 

Jo 

There  has  been  a  sustained  interest  in  proving  the  Central  Limit  Theorem,  for  t  — »■  oo, 
for  such  functionals,  that  is,  determining  the  limiting  distribution  of  the  normed  random 
variable 

(12)  j‘f(X(a))ds-EUinX(,))<U) 

{Far  f(X(s))  ds yn 
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A  CENTRAL  LIMIT  THEOREM  FOR  THE  RENORMALIZED 
SELF-INTERSECTION  LOCAL  TIME  OF  A  STATIONARY 

VECTOR  GAUSSIAN  PROCESS1 

By  Simeon  M.  Berman 
New  York  University 

Let  X(t)  be  a  stationary  vector  Gaussian  process  in  R"1  whose  compo¬ 
nents  are  independent  copies  of  a  real  stationary  Gaussian  process  with 
covariance  function  r(t).  Let  4>(z)  be  the  standard  normal  density  and,  for 
t  >  0,  e  >  0,  consider  the  double  integral 

m 

nte-m 

0  0  j»l 

which  represents  an  approximate  self-intersection  local  time  of  X(s),  0  s  s 
£  t.  Under  the  sole  condition  r  e  L2,  the  double  integral  has,  upon  suitable 
normalization,  a  limiting  normal  distribution  under  a  class  of  limit  opera¬ 
tions  in  which  t  -» o°  and  e  =  e(i)  tends  to  0  sufficiently  slowly.  The 
expected  value  and  standard  deviation  of  the  double  integral,  which  are  the 
normalizing  functions,  are  asymptotically  equal  to  constant  multiples  of  t2 
and  t3/2,  respectively.  These  results  are  valid  without  any  restrictions  on 
the  behavior  of  r(t)  for  t  -»  0  other  than  continuity. 


1.  Introduction  and  summary.  Let  X(t),  t  >  0,  be  a  real,  measurable 
stationary  Gaussian  process.  For  simplicity,  take  EX(t)  =  0  and  EX2(t)  =  1 
and  let  r(t)  =  EX(0)X(t)  be  the  covariance  function,  which  is  assumed  to  be 
continuous.  For  m  ^  1,  let  XJit), : . . ,  Xm(t),  t  >  0,  be  independent  copies  of 
X(t),  and  define  the  vector  process  X(t)  =  (X^t), Xm(t)).  Put 


(1.1) 


and,  for  e  >  0  and  t  >  0,  consider  the  random  variable 


(1.2) 


Xj(s)  -  Xj(s') 


dsds'. 


The  following  theorem  is  our  main  result. 
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THE  TAIL  OF  THE  CONVOLUTION  OF  DENSITIES  AND  ITS 
APPLICATION  TO  A  MODEL  OF  HIV-LATENCY  TIME1 

By  Simeon  M.  Barman 
New  York  University 

Let  p(x )  and  q(x)  be  density  functions  and  let  (p  *  qXx)  be  their 
convolution.  Define 

w(x)  —  —  (d/dx)logq(x)  and  v(x)  =  —  (d/dx)log  p(x). 

Under  the  hypothesis  of  the  regular  oscillation  of  the  functions  w  and  u, 
the  asymptotic  form  of  (p  *  qXx),  for  x  -»  °°,  is  obtained.  The  results  are 
applied  to  a  model  previously  introduced  by  the  author  for  the  estimation  of 
the  distribution  of  MV  latency  time. 

1.  Introduction  and  summaiy.  Let  p(x)  and  q(x)  be  probability  den¬ 
sity  functions  and  (p  *  qXx)  their  convolution.  The  focus  of  this  paper  is  the 
determination  of  the  asymptotic  form  of  ( p  *  qXx)  for  x  -*  «  or  x  -*  6  =  sup 
(support  of  p  *  q)  on  the  basis  of  the  asymptotic  forms  of  p(x)  and  q(x).  In 
most  of  this  paper  p  and  q  are  assumed  to  be  of  different  orders  of  magnitude 
for  x  -*  oo  and  q  will  have  the  role  of  the  density  with  the  heavier  tail.  The  key 
tools  for  densities  p  and  q  with  unbounded  support  are  the  functions  v(x)  = 
—  (.d/dx)log  p(x)  and  w(x)  =  —(d/dx)\og  q(x),  which  are  closely  related  to 
the  hazard  functions  used  in  extreme  value  theory.  The  tail  of  q  obviously 
dominates  the  tail  of  p  whenever  the  reverse  holds  for  their  corresponding 
functions  w  and  v.  Throughout  this  paper  it  is  assumed  that  v(x)  and  w(x) 
are  nonnegative  for  all  sufficiently  large  x.  In  particular,  it  follows  that  the 
corresponding  densities  are  nonincreasing  for  such  x. 

Theorem  3.1  states  that  if  limsup  w(x)  <  liminf  v(x)  for  *-*<»,  then 

/  P(*  ~  t)q(t)  dt  ~  q(x)  f  etwWp(t)  dt. 

This  represents  an  extension  of  a  corresponding  result  of  Breiman  (1965)  and 
Cline  (1986),  stated  in  terms  of  distributions  instead  of  densities  for  the  case 
w(x )  ->  e  >  0.  Theorem  3.2  furnishes  a  general  result  under  the  condition 
w(x)/v(x)  -*  0,  which  allows  for  even  the  cases  w(x)  -*  oo  or  v(x)  ->  0.  It 
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A  CENTRAL  LIMIT  THEOREM  FOR  THE  RENORMALIZED 
SELF-INTERSECTION  LOCAL  TIME  OF  A  STATIONARY  PROCESS 


Simeon  M.  Berman* 


1.  Introduction  and  Summary 

Let  X(t),  t  >  0,  be  a  measurable,  separable  stochastic  process  in  Rm,  for  some 
m  >  1.  Let  4>(x),  x  €  Rm,  be  a  probability  density  function  such  that  4>{x)  =  <f>(—x), 
and  let  4>(u)  be  the  corresponding  characteristic  function.  It  is  assumed  that 

(1.1)  /  |<£(ii)|  du  <  oo  . 

jRm 

Let  B  be  a  symmetric  (B  =  — B )  closed  bounded  Borel  set  in  Rm.  For  e  >  0  and 
t  >  0  define  the  functional 


(1.2)  I.(t)  «  £  jf  «—  <S  (£W_£M)  lB(X(s))  lB(X(s'))dsds'  ; 


this  is  the  approximate  self-intersection  local  time  of  X(s),  0  <  s  <  t,  relative 
to  B,  for  small  e  >  0.  (The  ratio  X(s)  —  A(s')/e  in  (1.2)  is  understood  as  e-1 
times  the  vector  X(s)  —  X(s').)  This  is  a  more  general  version  of  the  functional 
introduced  by  the  author  (Berman,  1992)  in  the  context  of  a  stationary  Gaussian 
process  X(t)  in  the  special  case  where  B  =  Rm  and  <j)  is  the  m- dimensional  prod¬ 
uct  standard  normal  density.  The  inclusion  of  the  indicator  1b  in  the  integrand 
restricts  the  approximate  self- intersections  to  those  for  which  the  sample  function 
values  belong  to  B.  In  the  previous  work  it  was  shown  that  if  X(t ),  t  >  0,  is  a 
stationary'  Gaussian  process  in  Rm  whose  components  are  independent  copies  of  a 
real  stationary  Gaussian  process,  then,  under  a  mild  condition  on  the  covariance 
function,  the  random  variable  t~3t2(I,(t)  —  EIf{t) )  has,  under  an  appropriate  class 
*  Supported  by  NSF  grant  DMS  SS-0118S  and  ARO  contract  DAAL  03-S9-K- 
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of  limit  operations  t  — +  oo,  e  =  e(t)  — *  0,  a  limiting  N( 0,  a2)  distribution,  where  a2 
is  explicitly  given.  The  current  work  extends  this  central  limit  theorem  to  the  func¬ 
tional  (1.2)  and  where  the  class  of  stochastic  processes  X(t)  includes,  in  addition  to 
the  Gaussian  process  studied  in  Berman  (1992),  a  large  class  of  ergodic  stationary 
Markov  processes. 

Our  main  result  is 

THEOREM  1.1.  Let  X(i)  be  a  stationary  vector  process  in  Rm  with  marginal 
density  p(x),  x  €  Rm.  Let  p(x,y\t)  represent  the  joint  density  of  X{ 0)  and  X{t)  at 
(x,y),  for  t  >  0.  Let  t  —  e(t)  be  a  decreasing  function  for  t  >  0  with  e(t)  — *■  0  for 
t  — *  0,  where  the  convergence  is  so  slow  that 

(1.3)  kin  t~1/2(e(t))~m  J  J&  j  <t>  p(x ,  y;  s)  dx  dy  ds  =  0  . 

Assume: 

(1.4)  p(x)  is  continuous  on  B  , 
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for  t  —*  oo.  In  this  case  a2  is  given  by 

(1.10)  a2  =  2  f  [  p(x)p(y)  [  [p(x,y,s)-p(x)p(y)]dsdxdy  . 

Jb  Jb  Jo 

We  note  that  the  condition 

(1.11)  il/2em(t)->  oo 

is  sufficient  for  (1.3)  because  <f>  is  bounded  under  (1.1)  and  f  f  p(x,y;s)  dx  dy  =  1 
for  every  s  >  0.  We  also  note  that  the  scaling  factor  t3^2  in  (1.9)  is  the  same  for  all 
dimensions  m  >  1. 

For  c  — *  0  the  integral  (1.2)  measures  the  set  of  points  (s,  s')  where  A’(s)  and 
X(s')  are  close.  This  includes  both  points  near  the  diagonal  and  points  bounded 
away  from  it.  The  latter  are  the  “genuine”  near-self-intersection  points.  We  show 
that  the  near-diagonal  points  make  an  asymptotically  negligible  contribution  to 
I((t)  in  the  statement  of  Theorem  4.1  in  the  sense  that  (1.9)  still  holds  after  the 
removal  of  these  points  from  the  domain  of  integration  in  (1.2).  For  this  purpose, 
we  prove: 

THEOREM  1.2.  Under  ihe  conditions  of  Theorem  1.1, 


(1.12) 


,LV3/’£(//  1[o,i]<l*-''l> 


^(£i£L_£if!i)lB(x(5))iB(AV))*d5')  =  0 . 


PROOF.  The  expected  value  in  (1.12)  is  at  most  equal  to  (see  Section  6) 


2J  J  1[0,1  ~  jB<t>{—^-)p(x,y-,s' -  s)dxdyds' ds 

(  <2 1  [  e~m  [  [  <t>(- — -)p(x,y;s)dxdy  ds  , 

Jo  Jb  Jb  € 

which,  under  the  condition  (1.3)  is  of  order  smaller  than  t 3^2.  □ 

The  application  of  Theorem  1.1  to  Markov  processes  X(t)  is  demonstrated  in 
Section  5.  The  bivariate  density  p{x,  y;  s)  assumes  the  particular  form  p(x )  q(s;  x,  y), 
where  q  is  the  transition  density.  The  hypotheses  (1.4)-(1.7)  can  be  stated  in  terms 
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ABSTRACT 

LetX(r),r  >  0,  be  a  real-valued  diffusion  process  having  a  stationary  probability  measure. 
For  an  increasing  function  u(t),t>  0,  put  L(t)  =  mes  (s:0<s  <t,X(s)  >  u(t)X  It  is 
shown,  under  general  conditions  on  the  diffusion  coefficients,  that  if  u(t)  «  at  a 
sufficiently  slow  rate,  then  (L(f)  -  EL(t))l(VarL(t)),/l  has,  for  t  — » «,  a  limiting  normal 
distribution.  The  rate  of  increase  of  u(t)  is  stated  in  terms  of  the  scale  function  S(x)  asso¬ 
ciated  with  the  generator  of  the  process;  u(t)  must  satisfy  =  o(r),  for  t  ->  °°.  This 

complements  an  earlier  result  (Berman,  1988)  in  the  case  S(u(t ))  ~  r,  where  it  was  shown 
that  there  is  a  function  v(r)  such  that  v(r)L(r)  has  a  particular  limiting  compound  Poisson 
distribution. 
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